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ABSTRACT 

We investigate the variability of C IV A1549 broad absorption line (BAL) 
troughs over rest-frame time scales of up to ~7 yr in 14 quasars at redshifts 
z > 2.1. For 9 sources at sufficiently high redshift, we also compare C IV and 
Si IV A 1400 absorption variation. We compare shorter- and longer-term vari- 
ability using spectra from up to four different epochs per source and find com- 
plex patterns of variation in the sample overall. The scatter in the change of 
absorption equivalent width (EW), AEW, increases with the time between ob- 
servations. BALs do not, in general, strengthen or weaken monotonically, and 
variation observed over shorter (<months) time scales is not predictive of multi- 
year variation. We find no evidence for asymmetry in the distribution of AEW 
that would indicate that BALs form and decay on different time scales, and we 
constrain the typical BAL lifetime to be >30 yr. The BAL absorption for one 
source, LBQS 0022-1-0150, has weakened and may now be classified as a mini- 
BAL. Another source, 1235-1-1453, shows evidence of variable, blue continuum 
emission that is relatively unabsorbed by the BAL outfiow. C IV and Si IV 
BAL shape changes are related in at least some sources. Given their high ve- 
locities, BAL outfiows apparently traverse large spatial regions and may interact 
with parsec-scale structures such as an obscuring torus. Assuming BAL outfiows 
are launched from a rotating accretion disk, notable azimuthal symmetry is re- 
quired in the outfiow to explain the relatively small changes observed in velocity 
structure over times up to 7 yr. 
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Subject headings: galaxies: active — galaxies: nuclei — X-rays: general — 
quasars: absorption lines — quasars: emission lines — quasars: individual (LBQS 
0009+0219, LBQS 0019+0107, LBQS 0021-0213, LBQS 0022+0150, LBQS 0025- 
0151, LBQS 0029+0017, LBQS 1231+1320, LBQS 1235+0857, LBQS 1235+1453, 
LBQS 1240+1607, LBQS 1243+0121, LBQS 1314+0116, LBQS 1442-0011, LBQS 
1443+0141) 



1. INTRODUCTION 



High-velocity, structured outflows in the central regions of quasars (QSOs) generate 
broad absorption lines (BALs) from material along the line of sight to the QSO ultraviolet 
(UV) emitter. These "BAL outflows" are launched, ionized, and accelerated by physical 
processes in the heart of the QSO, and can provide an important means of carrying ma- 
terial and energy out of the QSO's central region. BAL features can reach high outflow 
velocities (30,000 km s"^ or more), and t he absorption signatu res of these outflows are, 
by definition, at least 2000 km s~^ wide (jWeymann et al.lll99l[ ). although BALs can be 
much broa der. BALs are observed in up to ~15% of QSOs in optical/UV spectroscopic sur- 



veys 


(e.g., 


Wevmann et a 


2006 




Gibson et al. 


2009b 



19911 : iHewett fc Foltzl 120031 : iReichard et all l2003l : iTrump et al. 



and references therein), and the fraction of QSOs with BALs 



rises to 17-23% after correcting for optic al/UV selection effects (e.g.. iHewett fc Foltzl 12003 



Knigge et al.l l2008l : iGibson et al.l l2009bl ). The intrinsi c fraction may be ^2 times higher 



if less-restrictive criteria are used to select BALs (e.g., iTrump et al.ll2006t iDai et al.l 12008 



Ganguly fc BrothertonI 120081 ) ; in particular, narrower, mini-BAL absorption features ma y 
also be related to BALs (e.g., lArav et al.lll999b[ lOallagher et al.ll2002[ loibson et al.ll2009al ). 
Some physical models associat e BAL outflows with equatorial accretio n disk winds that may 
be ubiquitous in QSOs (e.g., iMurray et al.l Il995t iProga et al.l l2000l ). BAL outflows may 

vhen QSOs expel thick shrouds of gas and dust (e.g.. 



Voit et al. 


1993; 


Greee et al. 


2006) 



High-resolution spectroscopic studies have found that BAL absorption i s frequ ently 
saturated, but does not completely obscure the continuum emitter. lArav et al.l ( 1l999al ) con- 
cluded that BAL absorption troughs are generally shaped by the outflow kinematic structure 
and geometry, with the depth of the absorption feature representing the fraction of the con- 
tinuum emitter covered by one or more components of the absorbing outflow at a given 
velocity. In this model, studies of BAL variability primarily reveal ongoing changes in the 
structure of the outflow. Variability studies have also sought to test the "covering factor" 
hypothesis by identifying correlations between absorption and emission variation that would 
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indicate that absorption is not strongly saturated, but that the absorption strength is related 
to the photoionization-dependent column density of the absorbing ion. The results of these 
studies are mixed, but conclusive evidence for pho toionization-dominated variation has not 
been found in the general p opulation of BALs (e.g 
2OO7I : ICibson et alJboOSbf ). 



Barlow et al.lll989l . Il992l : iLundgren et al. 



With sufficient spectroscopic observing resources, BAL monitoring studies could greatly 
improve our knowledge of the composition and evolving structure of QSO outflows, and 
the physical mechanisms that launch and accelerate outflows. In practice, however, BAL 
variability studies have been limited by several constraints. Because BAL absorption is 
most commonly studied in the C IV A1549 line, BAL observations must be obtained either 
in the UV or of sources at sufficiently high redshift to place C IV BALs into a bandpass 
observable from ground-based telescopes. Because most known BAL QSOs are at z > 1.5, 
it can be difficult to study variation over extended (rest-fram e) time scales. W ith modern 
AGN surveys such as the Sloan Digital Sky Survey (SDSS; lYork et all I2OOOI). it is now 



possible to identify thousa nds of QSOs with observed BAL absorption (e.g., iTrump et al. 



20061 : iGibson et al.ll2009bl ). Yet it is necessary to have spectra obtained > 10 yr apart in 



order to span several rest-frame years. 

In this work, we examine the variation of BAL QSOs observed over rest-frame time scales 
of up to 7 yr. In order to construct a reasonably r epresentative sample, w e have selected all 
BAL QSOs at redshifts z > 2.1 fro m the catalog of Hewett fc Foltzl (120031 ) that are visible to 



the Hobby-Eberly Telescope (HET; iRamsey et al 



19981 ). Our s ample therefore inh erits the 



selection properties of the Large Bright Quasar Survey (LBQS; iHewett et al.lll995l ). We re- 
observed each of the 14 QSOs match ing our selection criteria using the HET Low Resolution 
Spectrograph (LRS; iHill et al.lll998l ) in order to obtain spectra with adequate signal-to-noise 
ratio (S/N) and spectral resolution for comparison to the LBQS epoch. Additionally, many 
of the sources in our study have s pectra available f rom i ntermediate epochs, obtained either 
as part of the follow-up study by lWeymann et al.l (Il99l[ ). or during the course of the SDSS. 
We will refer to the former set of spectra as the "Palomar" BAL spectra. The LBQS spectra 
in our current study and that of iGibson et al.l (l2008bl . hereafter G08) were obtained in the 
years 1986-1989. With the addition of the Palomar and SDSS spectra, we can sample time 
scales of months to several years for individual sources. 

In general, BALs vary in complex ways over multi-year time scales (G08). While studies 



of individual BAL QSOs are needed to und erstand how BALs appe ar (e.g., iHamann et al. 



of individual BAL ( 


20081: 


Leiehlv et al. 


Hall et al. 


2007 


).it 



Leighlv et al.l l2009l ). disappear (e.g.. lJunkkarinen et al.l l200ll ). and accelerate (e.g. 



20071 ). it is also important to determine the general characteristics of BAL outflow 



evolution. In the current work, we take another step in that direction by comparing C IV 
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BAL absorption profiles in up to four epoclis spanning a range of time scales from months to 
7 yr. Many of our sources are at sufficiently high redshift to shift the Si IV absorption region 
into the optical bandpass, enabling a study of the simultaneous variation of absorption from 
two different ions (Si IV and C IV). 



1.1. Conventions 

Because conventions vary among studies, we briefly describe the terminology used in 
this work. We define outfiow velocities to be negative (with respect to QSO emission rest 
frames). Positive velocities would indicate features that are at longer wavelengths than 
the wavelength corresponding to (rest-frame) zero velocity. However, we use the terms 
"greater" and "smaller" velocities to refer to the magnitude of the velocity, so that an outfiow 
speed of -10,000 km s~^ is "greater" than a speed of -5000 km s~^. We consider absorption 
equivalent widths (EWs) to be negative, so that a negative change in AEW between two 
epochs corresponds to a strengthening BAL. Errors in EW are calculated by propagating 
standard errors, with an estimate of 10% error in the continuum. 

Although we often quote the averaged doublet wavelengths to identify a spectral line, 
we perform absorption outfiow calculations using the wavelength of the red component of the 
doublet. In effect, we use absorption rest wavelengths of 1550.77 A for C IV and 1402.77 A 
for Si IV. This convention was adopted because the onset of a strong absorption trough from 
outfiowing material would be associated with the red side of the doublet. Unless otherwise 
noted, wavelengths in this work refer to rest-frame values. We use vacuum wavelengths for 
all epochs. 

Sources identifiers that are formatted like "0009-1-0219" refer to LBQS source des- 
ignations using B1950.0 coordinates. Throughout, we use a cosmology in which Hq = 
70 km s"^ Mpc-\ Qm = 0.3, and = 0.7. 



2. OBSERVATIONS AND DATA REDUCTION 



2.1. Observations and Methods 



Our sample is c ompr ised of the 14 BAL QSOs at 2; > 2.1 identified in the LBQS survey 
by iHewett fc Foltzl ( l2003l ) with C IV absorption regio ns that can be observed by the HET. 
We selected sources from the iHewett fc Foltzl (|2003[ ) list with redshifts 2; > 2.1 so that 
potential high-velocity absorption extending blueward to 1400 A could be observed in the 
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HET LRS bandpass (observed-frame wavelengths >4340 A). One source, LBQS 2212-1759, 
was excluded because its declination was lower than the HET limit. The observations used 
in our study are listed in Table [T] 

LBQS spectra were obtained from 1986-1988 with the Multiple Mirror Telescope (MMT) 
spectrograph. They cover a wavelength range of 3200-7400 A and g enerally have a sp ectral 
resolution of ~6 A with wavelength calibrations accurate to ~4 A ( Foltz et al.l 119871 ). The 
"Palomar" spectra we re taken as part of a campaign to obtain spectra of LBQS QSOs, 
especially BAL QSOs ( IWeymann et al.lll99ll ). Most of the Palomar spect ra were obtained 
using the Hale Telescope with the Double Spectrograph ( lOke fc Gunnlll982[ ) from 1988-1989 
with a resolution of about 6.5 A in bluer wavelengths. These spectra are available for most of 
our sources. We discarded three cases (0021-0213, 0025-0151, and 0029+0017) for which the 
original spectrum was not available at rest-frame wavelengths A < 1550 A. We also omitted 
the spectrum for 1231+1320 because of apparent inconsistencies in its flux correction at 
bluer wavelengths. SDSS sp ectra for four sources were obtai ned from the public SDSS Data 
Release 5 (DR5) database (lAdelman-McCarthy et al.l 120071 ). These spectra have a higher 
resolution (~3 A) than the other spectra in our study. 

Since 2007, we have also been observing these QSOs with the HET LRS. We used the 
g2 grism and a GG385 filter with an 1.5" slit, giving wavelength coverage down to ^4340 A 
on the blue side and a spectral resolution (estimated from fits to sky lines) of ^6.2 A, similar 
to those of the LBQS and Palomar spectra. We reduced the spectra using standard IRAF 
procedures and fiux-corrected each spectrum using a standard star observation obtained the 
same night. One of our sources, 1235+1453, had already been observed by our group in 
2002, so we did not re-observe that source. Exposure times for our sources ranged from 
~ 15-20 minutes. In a small number of cases, we interpolated over bins to remove "glitches" 
in the spectra that were much narrower than the instrumental resolution and are typically 
associated with cosmic ray hits or detector defects. 

Throughout this work, we focus on absorption (and emission) profiles obtained by divid- 
ing observed spectral features by an estimated model of continuum and broad-line emission. 
Given the possible errors in flux calibration for the HET and other observing epochs that 
are beyond our control, we do not depend on absolute or relative flux calibration. Addi- 
tionally, visual inspection indicated that the wavelength calibration between telescopes may 
have differed slightly for up to 5 sources. We examined non-BAL absorption and emission 
in spectral regions surrounding the broad C IV emission line at 1549 A to visually deter- 
mine a multiplicative correction factor (<0.3% in 4 of 5 cases) that we applied to LBQS 
wavelengths to obtain a satisfactory match among epochs. We selected a fiducial common 
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redshift for every epoch of each source from the SIMBAD databasejl] these values are hsted in 
Table [H These small calibration adjustments are sufficient for our study, which relies upon 
a relatively small wavelength range (1300-1550 A) and does not require precise absolute 
wavelength calibration. 

For each spectrurn in our study, we fit an emission model following the prescription in 



G08 and iGibson et al.l ( l2009bl ). In brief, we fit an intrinsically reddened power law contin- 
uum model to each spectrum after correcting for Galactic reddening. We initially fit only 
wavelength regions that are expected to be relatively free of emission or absorption features. 
We then iteratively fit emission lines to Si IV A1400, C IV A1549, Al III A1857, C III] A1909, 
and Mg II A2799. At each stage, we exclude regions that deviate strongly from the fit in 
order to exclude strong absorption and emission features. We model the cores and wings 
of each emission line using a single Voigt profile for convenience, although we attach no 
physical significance to this model. Our approach to fitting emission line regions is to fit 
the emission that is present, excepting regions that exhibit obvious absorption features. If 
emission lines are absorbed by variable absorbing components that are broad and smooth 
enough to cover the broad emission line without leaving obvious absorption features, we 
would not be sensitive to this effect. Of course, emission lines may also vary for a number 
of reasons unrelated to absorption. Also, any variable asymmetry in broad emission line 
profiles will not be perfectly modeled with our symmetric line profiles. Our chosen approach 
is motivated partly by the capabilities of the current data set; future data sets with improved 
fiux calibration and spectral resolution will be able to address these issues in more detail. 

We inspected and manually adjusted the fits in each case, using spectra from other 
epochs as a guide in a small number of cases where there was ambiguity in continuum 
placement or emission line shape. When adjusting continua, we ignored the region from 
~2200-3000 A that can be contaminated by ionized Fe emission; such emission is difficult to 
constrain and the continuum at those wavelengths does not strongly infiuence our estimates 
of Si IV and C IV properties in any case. The spectra and emission models for each source 
and epoch are shown in Figure [H 

Most of the data a nalysis for this project was performed using the ISIS platform 



f lHouck fc Denicolall2000l ) 



f 



http : //simbad . u-strasbg . f r/simbad/sim-f id 
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2.2. The Case of 1235+1453 



One of the sources, 1235+1453, is an outher in several of our analyses. The C IV BAL 
of this source is dramatically weaker in the Palomar epoch. Figure [2] compares the LBQS, 
SDSS, and HET spectra to the Palomar spectrum of 1235+1453. The figure shows that the 
continuum of 1235+1453 has also changed dramatically over time; it is much bluer in the 
earlier (LSST and Palomar) epochs. This is most clear in the middle panel of Figure |2l 
where we compare the two epochs with the best fiux calibration (Palomar and SDSS). 

The values of the C IV BAL EW at the different epochs depend on the continuum level 
in each spectrum. We have experimented with different continuum fits to model the Palomar- 
epoch BAL absorption and are unable to find a satisfactory model that matches the apparent 
blue continuum and also results in a C IV BAL resembling those in other epochs. Figure [3] 
demonstrates two alternative continuum placements. For this work, we have opted to use 
the higher continuum placement. The lower continuum attempts to match the continuum 
level "underneath" the emission line at 1300 A, but also falls below the emission levels at 
1400-1450 A and in the absorbed region blueward of 1200 A. An even higher continuum 
placement would be required to obtain a BAL strength equivalent to other epochs. 

The existing photometry and spectra for 1235+1453 are consistent with the presence 
of a very blue ultraviolet continuum component that declines over time and is not absorbed 
by the BAL outfiow. The differences in the Palomar-epoch BAL are due to changes in 
the absorption trough depth (by 20-25% of the continuum level), rather than changes in the 
absorption profile. A variable, blue continuum component would also explain the substantial 
changes in the (rest-frame) UV spectral shape of 1235+1453 as well as the fainter gf-band 
magnitude at the epoch of the SDSS spectrum relative to that of the Palomar spectrum 
(^0.6 mag fainter over a rest-frame interval of ^4 yr). From the current data, the cause of 
the observed blue variability is not obvious. SD SS spectroscopic st udies have demonstrated 
that quasars are generally bluer when brighter (jWilhite et al.ll2005r). and the spect ral shape 
variability could be associated with changes in accretion rate (jPereyra et al.ll2006l ). 



The amplitude of the variability cha nge shown is unusual bu t not extreme for a 



repre- 
sentative quasar over several years (e.g., IVanden Berk et al .112 0041 : IWilhite et al.l 2005). W e 



also note that 1235+1453 was also somewhat unusual in the study of I Gallagher et al.l ( l2006l ). 
in that it showed the maximum observed degree of X-ray weakness out of a sample of 35 
LBQ S BAL QSOs. While the variable continuum of quasar spectra has been studied in detail 



(e.g.. iWilhite et al.ll2005l ). we believe our discovery that the variable continuum component 
of 1235+1453 does not show the presence of high-ionization BAL troughs is a new result. 
Confirmation of such behavior among the BAL QSO population would provide significant 
new constraints on the origin of BAL outfiows and their geometry with respect to quasar 
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emission regions. 

Although the behavior of 1235+1453 is unusual in the context of the BAL QSO sample 
investigated here, we do not wish to bias our data by removing 1235+1453 from our sample. 
We retain 1235+1453 in our study, but we note cases where it is an outlier that significantly 
affects our results. 



3. ANALYSIS AND DISCUSSION 

FigurelHand FigureEJshow the "ratio spectra" of the Si IV and C IV absorption regions, 
respectively, with our emission model divided out. The figures demonstrate that patterns 
of BAL variation are, in general, complex. Rather than attempting to describe the wide 
variety of changes in these absorption features, we focus instead on statistical statements 
about the evolution of BAL absorption. We rely primarily on measures such as the EW that 
are integrated over the whole wavelength range where absorption could occur. We prefer 
such measures to thos e that are calculated o n a set of bins specific to each BAL, such as the 



"balnicity index" (BL iWeymann et al.lll99ll ) and maximum outflow velocity (fmax), because 
these measures are highly dependent on the wavelengths chosen to be the boundaries of 
each BAL trough. BAL boundaries can change dramatically due to spectral resolution and 
random fluctuations in even a single spectral bin, and we do not wish the results of our study 
to be unduly influenced by our choice of BAL boundaries in each observing epoch. 

Of course, the integrated measure of EW is not sensitive to changes in which one part of 
a BAL weakens, while another part strengthens. In some parts of our study (e.g., §3.51 and 
we investigate the patterns of variation in more detail. Additional statistics could also 



be employed to measure, e.g. , the sum of the magnitude of change in absorption profiles. Such 
statistics may be particularly useful in larger-scale studies using spectra having uniformly 
high resolution and S/N. 

Throughout this work, we search for correlations between absorption and emission prop- 
erties using the Spearman Rank Correlation test. In addition to describing each test in the 
work below, we summarize the results of our correlation tests in Table [2l We do not generally 
draw conclusions from correlations at a confidence level < 99%. 



3.1. Variation in C IV Equivalent Width 



We calculate the rest-frame EW of each absorption region using the continuum and 
emission-line fits described in §2.11 To calculate C IV absorption EWs, we integrate over 
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the entire region blueward of the C IV rest wavelength (1550.77 A) out to a velocity of 
-30,000 km s~^. We integrate EWs in units of Angstroms or km s~^; a conversion of 
^200 km s~^ per Angstrom is a good approximation for the C IV absorption region. Our 
calculated values are listed in Table Ei 

Figure [6] shows the evolution of EWs over time. For each source, we subtracted the ob- 
serving date and measured EW for the first (LBQS) epoch from those of the following epochs 
so that the LBQS BALs correspond to the point (0, 0) for each source. All measurements 
for a single source are connected with a solid line in the figure, and individual measurements 
are labeled according to the observing campaign ("P" for Palomar, "S" for SDSS, and "H" 
for HEX). 

In order to investigate whether the spread of AEW increases with time, we calculated 
s, defined as the square root of the unbiased sample variance of AEW, for the 8 sources 
with Atsys < 500 days and also for the 13 sources with Atsys > 2000 days. Here, Atgyg 
represents the elapsed time in the rest frame. For the first set, comprised solely of Palomar 
observations, s{Atsys < 500) ~ 1380 ±370 km s~^ and the median time between observations 
is 96 d. For the latter set, s{Atsys > 2000) ~ 1830±370 km s~^ and the median time between 
observations is 2 150 d. Here we estim ate the variance on as war(s^) ~ 2s^/ {N — 1) for 



data points (e.g.. lEdelson et al.ll2002 ): from this we estimate the error on s as s/ ^/2{N — 1] 



Our measured trend of increasing spread in AEW over time is qualitatively consistent 
with that of G08, although the difference in variance between the LBQS-Palomar and LBQS- 
HET observations is not highly significant and the variance for the Palomar set is somewhat 
larger than that measured over similar time scales in G08. The variance of the Palomar set is 
greatly increased by the single data point corresponding to 1235+1453 (see §2.11) . for which 
the BAL is dramatically weaker in the Palomar epoch. If we omit that outlying data point, 
the scatter of the set with Atgyg < 500 days decreases significantly to cr ~ 840 ± 240 km s~^. 
By contrast, removing the largest "outlier" from the set with Atgys > 2000 days decreases 
the variance by a much smaller amount, to 1610 ± 340 km s~^. We also investigated the use 
of the median absolute deviation to estimate the standard deviation a as: 

a ^ median(|xj — median(xj)|)/0.6745. 

The m edian absolute deviat ion is less sensitive to outlying points than the sample variance 



IS (e.g., 



Maronna et al.ll2006l ). In fact, we find a{Atsys < 500 d) = 880 km s ^ and cr{Atsys > 
2000 d) = 2100 km s'^, supporting the case that 1235+1453 is an outlier to the trend of 
increasing spread in AEW. Even so, it is not clear that 1235+1453 is appropriately classified 
as an outlier; much larger samples will determine the distribution of AEW and its apparent 
time evolution. 
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It may be the case that the distribution of AEW is not Gaussian over times spanning 
months to years. For example, sources may be moved to the tail of the AEW distribution 
by strong, temporary "outbursts" in continuum emission regions that are not covered by the 
absorbing outflow. Such outbursts could be independent of the outflow evolution. While 
the data for 1235+1453 suggest that the distribution of AEW may be governed by physical 
processes that make it non-Gaussian, our current sample is not large enough to verify this 
hypothesis. 

To analyze the general properties of long-term BAL variability over 3-7 rest-frame years, 
we combine the sample of this work with the (non- redundant) sources from G08. In each 
case, we compute AEW, the change in BAL EW between the earliest (LBQS) and latest 
(either SDSS or HET) epoch. A histogram of AEW values is plotted in Figure [71 For 
comparison, we have also plotted a histogram of AEW on the shorter time scales between 
LBQS and Palomar epochs for 8 sources using a dotted line in the figure. The distribution 
of AEW is reasonably symmetric, with 13 of 23 BALs having AEW < 0, corresponding to 
strengthening absorption. While there are subtle differences in source selection between the 
current work and GOB, and a weak bias toward strengthening BALs could be introduced 
by flux-based quasar detection in the earlier epoch, the overall result indicates no evidence 
of asymmetry in the increase or decrease of absorption EWs on multi-year time scales. 
Asymmetry could occur, for example, if BALs formed very rapidly and then decayed slowly; 
nearly all identified BALs would then be observed to decay over time. Any monotonic trends 
in BAL variation over 3-7 years must either be weaker than we can measure, and/or these 
trends are dominated by stronger variation on multi-year time scales. 



3.2. Variation as a Function of Absorption Strength 

In this section, we examine whether BAL variation depends on the overall strength (EW) 
of BAL absorption. There are several ways to do this. First, we look for relations between 
the amount of change in the C IV BAL (AEW or lAE'W^I) and a measure of BAL strength 
(LBQS-epoch EW or (EW), defined as the average EW between LBQS and HET epochs). 
We find no significant correlations (see Table [2]) excepting a relation (99.6% confidence) 
between lAE'l^l and the LBQS-epoch EW, and a weaker relation (97% confidence) between 
|Aii^iy| and (EW). These relations are not surprising, as stronger BALs are generally 
broader and have more absorbed regions that can vary. 



The amount of time between EW measurements may also affect the degree of change in 
EW. Following earlier studies (e.g.. iLundgren et al.ll2007l ). we define the fractional change in 
EW, \AEW/ {EW)\, where the average EW, (EW), is calculated by averaging the EW for 
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a given epoch with that of the LBQS epoch. (Note that we are technically discussing the 
magnitude of fractional variation, as we define this quantity to be positive.) In Figure IH we 
plot the magnitude of fractional change in EW as a function of the time interval between 
observations. We plot fractional changes in EW for the HET epoch and for the Palomar 
epoch, when available. For sources with Palomar spectra, we draw a solid line connecting 
the Palomar and HET epochs. As for the measurements of AEW, 1235+1453 is once again 
an outlier; it is the only source with a large decrease in \AEW/ {EW)\ between epochs. 



In agreement with previous studies (e.g., iLundgren et al.l 120071 : G08), we find that the 
largest fractional changes are observed on longer time scales. The cases with the largest 
fractional change in EW (0009+0219, 0025-0151, and 1314+0116) show reasonably strong 
BALs with BIq > 1000 km s~^ in one or both of the LB QS and HET epochs, where BIq is 
the traditional "balnicity index" ( iWeymann et al.l 1 199 if ) modified to include absorption in 
the region between —3000 and km s~^. 

As previously noted (G08), the large amount of variation observed in some strong BALs 
could, if continued at the same rate, lead to the disappearance of some BALs over ~15 
years. Because strong BALs seem to have lifetimes longer than that ( §3.41) . we would expect 
that the fractional EW variation "turns over" on longer time scales for these highly-variable 
sources and that they do not vary monotonically. We also note that 2 of the 3 sources with 
\AEW/ {EW)\ > 0.5 are actually increasing in strength, rather than decreasing. 



3.3. Variation on Sub- Year and Multi-Year Time Scales 

We also test for indications that changes in a significant number of individual BAL 
EWs may be monotonic over months to years. To do this, we compute AEWp^, defined 
as the difference in EW between the Palomar and LBQS epochs, and AEWhp, defined 
as the difference between HET and Palomar epochs. The typical rest-frame time between 
the LBQS and Palomar epochs is a few weeks to months, while the time between Palomar 
and HET epochs is 3.8-6.5 yr. We find no evidence of a correlation between AEWp^ and 
AEWi^p for the 8 sources in our sample that have spectra in all three epochs. While some 
BALs may vary monotonically over long time scales (including sources such as 1314+0116 in 
our sample), the variations observed in a BAL over a few months do not allow us to reliably 
predict the longer-term variation over several years for BALs in general. 

The rate of change in EW, AEW/ At, is also unrelated on shorter (LBQS to Palomar 
epochs) and longer (Palomar to HET epochs) time scales. We find no significant correlation 
in AEW/ At measured for these two time scales. In fact, several BALs show opposite signs 
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for A£^W^/At, as the absorption weakened between one set of observations and strengthened 
in the other (Figure |9]). 



3.4. Evolution of BAL Absorption 

In G08, we placed a loose lower limit on the lifetime of strong BALs of I^al > 18 yr? 
based on the observation that the strong BALs in the G08 study had not disappeared over 
typical rest-frame times of ~4.3 yr. G08 defined "strong BALs" to have Bio > 100 km s~^. 
All but one of the BALs in our current study have Big > for both the LBQS and HET 
epochs. The exception is 0022+0150, for which the absorption has weakened and BIq is 
measured to be zero in the HET epoch. In fact, the BIq measurement is extremely sensitive 
to continuum placement, and a small adjustment to the continuum would give a measurement 
of Bio > km s-\ 

We model the probabihty of observing a BAL to disappear to be Atgys/tBAL, where 
Atsys is the rest-frame time between observing epochs and tsAL is the typical BAL lifetime. 
Given that we have 19 strong BALs (with BIq > 100 km s~^ in at least one epoch) in our 
sample (including 5 non-redundant sources from GOB) and that the median Atgys is 5.8 yr 
for this augmented sample, we now find that we would have a >90% chance of observing 
more than one case (0022+0150) in which a BAL disappears if tsAL ^ 30 yr in the rest 
frame. (The lifetime rises to >50 yr if we allow 0022+0150 to be classified as a BAL in the 
HET epoch.) 

While BALs may be transient in individual cases, and weaker BALs (and mini-BALs) 
are not included in this estimate, it seems that strong BALs are overall rather stable on 
human time scales. BAL outflows therefore appear to travel large distances; material moving 
constantly at 10,000 km s~^ will traverse ~0.3-0.5 pc in 30-50 yr. Even if we are seeing 
a "standing outfl ow pattern" in which new material continually enters at an angle to our 



line of sight (e.g., Ide Kooll 119971 ). the material that has left our line of sight can still reach 
large distances over the BAL trough lifetime unless it is decelerated or obstructed after it 
has moved out of view. If BAL outflows are launched at accretion disk radii much smaller 
than 1 pc, they may still interact with any parsec-scale structures such as an obscuring torus 
(if present in luminous quasars at outflow inclination angles). 

The absorption for 0022+0150 has weakened to the point where the absorption trough 
widths are becoming narrower than the 2000 km s~^ cutoff required for BAL absorption. 
0022+0150 may be more properly classified as a "mini-BAL" QSO in the HET data and 



would be another demonstration of the link between BAL and mini-BAL QSOs (e.g.. lArav et al 
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1999bl : [Gallagher et al.ll2002l : iGibson et al.ll2009af ) . iLeighlv et al.l f l2009f ) have recently pointed 
out evolution in the opposite direction for the nearby Seyfert 1 galaxy WPVS 007 (at 
z ~ 0.03). In that source, BAL troughs emerged to augment the mini-BALs that were 
seen in an earlier epoch. 



3.5. Comparing Variation in Si IV and C IV BAL Shape 

Detailed comparisons of absorption profiles among i onization species ca n help to map 



the ionization structure of the BAL outflow. For example, lArav et al.l ( 1l999bl ) observed that 
the ionization-dependent structure of BAL components in PC 0946 + 301 was best explained 
if BAL outflows were structured to have lower-ionization "tubes" of material embedded in or 
associated with higher-ionization flows. The lower-ionization material could have a different 
(and smaller overall) covering factor than the higher-ionization material, leading to somewhat 
different trough shapes. As a consequence of the outflow geometry, lower-ionization material 
is also seen in a narrower velocity range. In fact. Si IV troughs do ap pear to be narrower 



than C IV troughs for BAL QSOs in general (e.g.. iGibson et al.ll2009bl ) 



Variability adds a new dimension to this picture. If absorption depths are determined by 
the fraction of the emitter covered by the outflow, then synchronous variation in absorption 
from different ions would require that the flow be structured so that different ionization 
components can vary in unison. On the other hand, if the absorption is not saturated, 
changes in the outflow ionization level may lead to more complex changes in the optical depth 
of the outflow. In this section, we discuss the relation between Si IV A 1400 and C IV A 1549 
BAL troughs in different epochs in order to determine whether they vary synchronously. Si IV 
and C IV have ionization potentials of 45.1 and 64.5 eV, respectively, and both ions have 
large ionization fractions over a range of (ionizing-continuum-dependent) absorber ionization 
levels. 

We use a straightforward approach in which we directly compare the observed changes 
in BAL absorption for Si IV and C IV. We do this using the four ratio spectra — LBQS 
and HET epochs for Si IV and C IV — that are obtained when the observed spectra are 
divided by our emission model. We smooth each spectrum using a boxcar filter of width 
3 bins and then rebin onto a common velocity grid. Our velocity grid only extends out to 
—20, 000 km s^^ to avoid contamination by unmodeled features blueward of Si IV. We then 
subtract the ratio spectrum of the earlier (LBQS) epoch from that of the later (HET) epoch. 
This gives us two "difference spectra," which we call AS'(f ) for the Si IV absorption region 
and AC(w) for C IV absorption. AS'(f ) represents the change in the Si IV absorption profile, 
and AC(w) does the same for C IV. Note that A5'(f) and AC(t') also include information 
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about the variation of narrower absorption features that may reside outside the BAL trough. 

Although the Si IV and C IV hne doublets have different widths that affect the different 
BAL trough shapes, this rough comparison still identifies correlations at >99.9% signifi- 
cance between AS'(w) and AC{v) for two cases, 0009+0219 and 1443+0141, according to 
a Spearman rank correlation test. Several other spectra show evidence for correlation at 
> 95% significance levels, although visual inspection indicates that the correlations for these 
may be influenced by a relatively small number of points on the boundaries of the AS{v) 
vs. AC{v) relation. At least 0009+0219 and 1443+0141 show coordinated changes between 
epochs over a range of velocities. Figure [TO] shows AS{v) and AC{v) for these sources. For 
a few other sources, the changes are suggestive or are strongly related over a smaller subset 
of the full velocity range. 



3.6. Comparing Variation in Si IV and C IV BAL Strength 

In this section, we are not concerned with the variation of absorption trough shapes, 
but only with the total strength of absorption. The absorption EW differs from ion to ion 
in a single epoch, w' ith Si IV absorption generally being weaker than C IV absorption (e.g.. 



Gibson et al.ll2009bl ). This could be due to a number of factors such as ionization levels, 
eleme ntal abundances, and BAL outflow geometry having an ion-dependent covering factor 
(e.g., 



Arav et al.lll999bf ). 



In Figure [TTl we plot AEW between LBQS and HET epochs for Si IV absorption 
regions compared to AEW for C IV absorption regions. Here again, we only integrate out 
to —20, 000 km s~^ for both ions in order to avoid emission model uncertainty for Si IV at 
higher outflow speeds. We only include the 9 sources at sufficiently high redshift to have 
HET coverage out to —20, 000 km s~^ for Si IV. Our calculated BAL EW values are listed 
in Table |H We find no significant correlation of AEW between Si IV and C IV. We also find 
no significant correlation of (EW) between the two ions. 

There is a suggestion of a correlation in the fractional change in EW, \AEW/ {EW)\, 
as shown in Figure [121 The correlation for the full set of sources is not highly significant, 
but a stronger correlation (> 99% confidence) would be observed if 1235+1453, the outlier 
at the top left of the plot, were removed. As previously discussed ( §2.ip . this source may 
have contaminating "emission" that artificially increases AEW and decreases (EW) , moving 
the data point for 1235+1453 upward in Figure [121 We use the IDL FITEXY code to fit 
F(SilV) against F(CIV), where F[x) is the fractional change in EW for ion x, accounting 
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for the formally- derived errors. We find a fit of: 

F(SilV) = (1.44 ±0.45)F(C IV) + (-0.05 ±0.19) (1) 

for the sample with 1235+1453 omitted. The fit is shown as a solid line in Figure O 

Correlations in \AEW/ {EW)\ can arise, for example, when changes in the ionization 
level or geometry affect the same fraction of material for each ion. While the ratios of Si IV 
to C IV may vary among different QSO outflows, it is plausible that the fraction of material 
that changes will be the same for both ions in a given outflow. For models of such scenarios, 
we would not necessarily expect to see a correlation between AEW for the two ions, as AEW 
can depend on the amount of absorbing material that is present along the line of sight for 
each ion (and also the precise mode of variation). 

Under the assumption that 1235+1453 can be classified as an outlier, our results would 
indicate that Si IV absorption varies in concert with C IV absorption, and in at least some 
cases, the changes in velocity-dependent absorption profiles can be highly related ( §3.5p . 
However, according to Figure [TT| it does appear that overall Si IV absorption can strengthen 
while C IV weakens, and vice versa. If variation were purely due to changes in the emission 
covering factor, we would expect that both ions would always strengthen or weaken together. 



3.7. Comparison With Emission Line Variation 



The strengt h (EW) of C I V line emission is known to be related to the UV continuum 
luminosity (e.g.. lBaldwinlll977l ). This "global Baldwin effect" relation apparently extends to 
higher-energy continuum emission as well, with the C IV emission EW being anti-correlated 
to the continuum mor iochromatic l uminosity at 2500 A and correlated to t he monochromatic 
luminosity at 2 keV JCreenl [iQQsl : IWang et allllQQsl : Icibson et aUboOSah . Additionally, an 
"intrinsic" Baldwin effect is observed in which C IV EW variation is anti-c orrelated with 
conti nuum luminosity in a single object as its spectrum varies over time (e.g., iKinney et al. 



1990l ). While BAL QSOs have weaker C IV line emission in general than non-BAL QSOs do, 
C IV emission in BA L QSOs still follows a "Baldwin effect" type trend with UV luminosity 
Jcibson et al.ll2009bh . 



Motivated by these considerations, we have briefly investigated whether we can use 
C IV emission EW as a proxy for the ionizing continuum luminosity to search for evidence of 
photoionization-driven BAL variation. We find no significant correlations between emission 
line EW (listed in Table [5]) and the C IV absorption EW. We also find no correlation between 
the changes in C IV emission EW and C IV BAL absorption. We do not draw strong physical 
conclusions from this result because any underlying correlation could be weakened by factors 
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such as different ionization time scales for tfie emitting and absorbing material as well as the 
fact that the observed C IV emission is integrated over a large spatial region. Furthermore, 
we are unable to account for any variable absorption that would have a velocity profile 
broad and smooth enough to affect the entire broad emission line without leaving obvious 
absorption signatures. 



3.8. The Stability of BAL Outflows 

While BAL shapes differ greatly from one QSO to another, numerous cases in the 
current study and also G08 demonstrate that the factors shaping BAL outflows (such as the 
velocity-dependent covering factor) remain relatively stable or even unchange d for years ove r 



much of the outflow velocity range. As previous researchers have noted (e.g.. lde Koollll997l ). 
BAL outflow models are challenged to explain both the range of variation effects seen in the 
spectra and also the remarkable stability in many velocity components of BAL outflows. In 
some cases, the velocity structure may be stabiliz ed by physical effects associated with the 



acceleration mechanism, such as line-locking (e.g.. iGallagher et al.ll2004l ) 



Over a median observed time of 5.8 yr in the rest frame, material flowing outward 
at 10,000 km s~^ travels a radial distance of 2 x 10^^ cm, which is 12,400 GMbh/c^ for 
a central black hole of mass Mbh = 10® ^0- If the BAL outflow is associated with a 
rotating accretion disk and resides at a radius of 10^^ cm, the Keplerian frequency of circular 
rotation is a; = \fGMBuT^ = 0.11 rad yr~^. This corresponds to a velocity of 3600 km s~^ 
perpendicular to the line of sight, or a distance of 7 x 10^^ cm covered over 5.8 years. Over 
a typical BAL lifetime of presumably >30-50 yr ( §3.4p . the azimuthal distances traveled 
would be >5-9 times larger, or > 0.1 — 0.2 pc. 

If the BAL outflow is launched from a rotating accretion disk, it must be possible to 
maintain relatively stable outflow structures over such large time and size scales. This is 
true whether we are observing a "standing outflow pattern" (in which the outflow material 
is continually replenished) or the same absorbing material over time. In the latter case, 
the outflow structure must also be stable as the absorber travels large radial distances. 
For a quasar with Mbh = 10® the broad emission line region radius (estimated from 



reverberation measurements of Balmer lines) is ~ 1.7 x 10^^ cm ( iKaspi et al.ll2000l ). and 
higher-velocity material in the BAL outflow has presumably traveled beyond this distance 
between the LBQS and HET epochs. 
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3.9. BAL Variation Patterns 

In this work, we have used absorption measures integrated over the entire absorption 
region (e.g., the absorption EW) to describe the variation of BAL troughs. In an earlier work 
(G08), we developed metrics that could be used to characterize the properties of variation 
in narrower velocity regions of a single BAL. We do not apply those metrics in the current 
study because they require SDSS-quality spectra with higher resolution and S/N in at least 
one epoch in order to be effective. However, visual inspection indicates that the variation 
in our sample is generally consistent with the G08 metrics. The variable regions of a BAL 
are typically narrower subsets of the full absorption trough, and BALs can vary in some 
absorption regions while remaining relatively constant in others (e.g., 0019+0107 at 1430- 
1480 A, Figure 5b; 0021-0213 at 1460-1480 A, Figure 5c). Variation occurs at both lower 
(e.g., 0025-0151 at 1510-1530 A, Figure 5e) and higher outflow velocities (e.g., 1314+0116 
at 1420-1480 A in HET, Figure 51; 1442-0011 at 1430-1450 A in HET, Figure 5m). We 
see variation in shallower (e.g., 0019+0107, Figure 5b) as well as deeper (e.g., 0021-0213, 
Figure 5c; 1231+1320 at 1520 A, Figure 5g) troughs, and the amount of change in the 
absorption depth is commonly < 30% of the continuum. 

Cross-calibration of the wavelengths between the LBQS and HET epochs reveals two 
potentially interesting cases of changes in velocity structure. In the first case, 0022+0150 
(Figure 5d), there are three relatively weak (by BAL standards) absorption features at ~1430, 
1490, and 1520 A. The strongest feature, at ~1520 A, has weakened over time. Additionally, 
there is an absorption feature on the blue side of the Lya line (Figured]). There is no clear 
way to line up the emission and absorption features between the LBQS and HET epochs so 
these features match. If we line up the spectra so that the red side of Lya emission matches, 
then the absorption features have shifted 5 — 10 A blueward. Alternatively, we can match up 
most of the broader and narrower absorption features, but the red wing of Lya has moved 
redward and the weak BAL absorption at 1430 A has moved blueward by about 5 A over 
time. Because the cross- calibration is anchored by absorption features within 130 A on both 
the blue and red sides of the 1430 A feature, any wavelength calibration issues would have 
to be highly nonlinear to account for the shift of the 1430 A feature. 

The second candidate for absorption acceleration is 1240+1607 (Figure 5j). While 
narrower absorption features match up well at ~1440, 1525, and 1535 A, the broader trough 
from ~ 1490 — 1510 A has shifted blueward by about 5 A over time. In this case, the 
absorption in the Palomar epoch from 1510 — 1540 A also follows a somewhat different 
pattern than the LBQS and HET epochs, suggesting that the absorption variation in this 
source may be particularly complex. 

These examples demonstrate issues that complicate a definitive search for BAL acceler- 
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ation. If the trough components have truly accelerated in these cases, the change in velocity 
is ~1000 km s"^ over 5.6 yr (0022+0150) and 6.0 yr (1240+1607) for these two sources. The 
corresponding acceleration is a factor of ^4-18 larger than th at observed in other sources 
(jVilkoviskij &: Irwiru l200ll : iRupke et al.ll2002l : iHall et al.l 120071 ). Of course, we cannot rule 
out the possibility that the troughs have simply changed shape rather than accelerating. If 
absorption deepens on the previously-unabsorbed blue side and weakens on the red side of 
a trough, the absorption feature will appear to accelerate. 



4. SUMMARY AND CONCLUSIONS 



Here, we briefly summarize results in our study. 



1. Multi-year BAL variability is complex; there is no clear pattern of C IV absorption 
variation over 5-7 yr in our sample of 14 BAL QSOs at 2; > 2.1. 

2. The scatter in AEW increases with time to ~1800 km s~^ over 5.8 yr. We find no 
evidence for asymmetry in the distribution of AEW about km s~^ that would indicate 
that BALs strengthen and weaken on dramatically different time scales. 

3. In one source, 1235+1453, we find evidence for a variable, blue continuum component 
that is relatively unabsorbed by the BAL outflow. 

4. Although individual BALs may strengthen or weaken monotonically, BALs do not 
generally evolve monotonically over 3-7 yr. Variation on multi-month time scales does 
not predict variation over years in an obvious way. 

5. Based on the fact that we saw at most 1 BAL (in 0022+0150) formally leaving our 
augmented sample of 19 strong BALs, we estimate that at least stronger BAL features 
have a lifetime >30 yr. The absorption in 0022+0150 did not fully disappear; instead, 
the absorption weakened to mini-BAL status. 

6. The evolution of substructures in Si IV and C IV BAL trough shapes appears to be 
related in at least two cases in our sample. 

7. While AEW is not strongly correlated over 5.8 yr for Si IV and C IV BAL trough 
evolution, there does appear to be a strong correlation in the fractional change in 
EW for Si IV and C IV troughs. This result depends on the removal of one outlier, 
1235+1453, which (as we discuss) has unusual properties in several respects. 
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8. Given their apparently long lifetimes and high velocities, BAL outflows can cover large 
radial and azimuthal distances relative to the accretion disk center. The stability of 
many BAL troughs (in whole or in specific velocity ranges) over years is an important 
consideration for models in which BAL outflows are associated with rotating accretion 
disks. 

9. We cross-calibrated the wavelengths of spectra between epochs using narrower absorp- 
tion and emission features and identifled two cases (0022+0150 and 1240+1607) of 
C IV BALs that showed evidence for possible changes in velocity structure of a mag- 
nitude ~1000 km s~^ over ^6 yr. Follow-up is needed for these and other candidate 
sources to determine the nature of possible absorber acceleration. 



The current work and p revious studies of BAL variability in QSO samples (e.g.. iBarlow 



19931 : iLundgren et al.l 120071 : G08) demonstrate that there is still much to learn about the 
structure, evolution, and influence of BAL outflows, as well as the physics that launches and 
shapes them. Improved samples of BAL spectra are clearly needed to examine the nature 
of BAL variation and acceleration, the ionization structure and evolution of BAL outflows, 
links with other absorption phenomena such as mini-BALs, and the long-term evolution of 
BAL outflows. 

Given the practical issues involved with observing UV spectra over long time frames, 
we suggest some approaches that would be particularly advantageous for the ongoing study 
of high-velocity QSO outflows. Firstly, larger samples are needed. The SDSS archive will 
provide an excellent baseline for studies of BALs at higher resolution and S/N, and the 
SDSS-III surveys will provide an additional epoch of flux-calibrated spectroscopy for many 
SDSS BAL QSOs. 

Secondly, intensive monitoring of a set of variable BALs, ideally at higher resolution, 
could resolve interesting questions about BAL structure and how BALs evolve over shorter 
time scales. Among other things, spectra from such a campaign could reveal substructures in 
the outflow that vary in concert (in the BAL of a single ion or from multiple ions), and would 
also enable a search for rapid photoionization-driven changes in the outflow. Candidates 
for BAL acceleration, such as 0022+0150, 1240+1607, and SDSS J024221.87 + 004912.6 



(jHall et al.ll2007l ) should also be re-observed over months to years to test whether acceleration 



has continued. 

Ideally, we would like to determine a "structure function for BAL variation" that would 
be a useful input to models of quasar outflows. The computation of structure functions and 
other detailed variability metrics will become possible as large-scale spectroscopic surveys 
such as SDSS-III re-observe the large samples of BAL quasars identifled in earlier SDSS 
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surveys. 
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Table 1. Observation Log 



LBQS Name 




z 


LBQS Date 


Palomar Date 


SDSS Date 


HET Date 


HET Exposure 
(min) 


Max Afsys 

(yr) 


0009 


+ 0219 


2, 


,642 


1987 


Sep 25 








2009 


Jan 1 


20.0 


5.8 


0019 


+ 0107 


2, 


.123 


1987 


Sep 25 


1988 


Aug 9 




2008 


Oct 22 


20.0 


6.7 


0021 


-0213 


2, 


.293 


1988 


Sep 9 








2007 


Dec 5 


20.0 


5.8 


0022 


+ 0150 


2, 


.770 


1987 


Sep 25 








2008 


Oct 21 


17.0 


5.6 


0025 


- 0151 


2, 


.076 


1987 


Nov 25 








2008 


Oct 1 


17.0 


6.8 


0029 


+ 0017 


2, 


.250 


1987 


Nov 25 






2001 Dec 19 


2008 


Oct 21 


17.0 


6.4 


1231 


+ 1320 


2, 


.380 


1986 


Dec 29 








2009 


Mar 3 


20.0 


6.6 


1235 


+ 0857 


2, 


.890 


1986 


Jun 6 


1988 


May 14 




2008 


Jun 29 


17.0 


5.7 


1235 


+ 1453 


2, 


.694 


1986 


Dec 29 


1988 


May 14 


2005 Mar 13 


2002 


Jun 14 


15.0 


4.9 


1240 


+ 1607 


2, 


.360 


1988 


May 10 


1989 


Mar 30 




2008 


Jun 28 


20.0 


6.0 


1243 


+ 0121 


2, 


.808 


1987 


Apr 2 


1989 


Mar 31 


2001 Jan 19 


2008 


Mar 1 


17.0 


5.5 


1314 


+ 0116 


2, 


.698 


1988 


Mar 16 


1988 


May 14 


2001 Mar 16 


2009 


Feb 25 


20.0 


5.7 


1442 


-0011 


2, 


.226 


1988 


Mar 16 


1988 


May 13 




2008 


Jun 28 


16.0 


6.3 


1443 


+ 0141 


2, 


.450 


1988 


Mar 16 


1988 


May 14 




2008 


Jun 28 


16.0 


5.9 



Table 2. Correlation Tests 



Comparison 


Reference 


Spearman 


Number of 




Section 


Confidence 


Data Pts 


AEVK vs. LBQS EW 


su 


69 


14 


|AEW| vs. LBQS EW 


332] 


99.62 


14 


AEW vs. (EW) 


332] 


1 


14 


|AEW| vs. (EW) 


332] 


97 


14 


AEWpi vs. AEWhp 


333] 


68 


8 


A£H//At (PL) vs. AEW/ At (HP) 


333] 


74 


8 


0009+0219: C\S{y) vs. AC(v) 


333] 


> 99.99 


166 


1443+0141: A5(d) vs. AC(v) 


3331 


> 99.99 


157 


(EW) for Si IV vs. C IV 


3331 


62 


9 


AEW for Si IV vs. C IV 


3331 


23 


9 


\l^EWH^EW)\ for Si IV vs. C IV 


3331 


76 


9 


ASM// (Ely) for Si IV vs. C IV 


3331 


99.35 


8 


(outlier 1235+1453 removed) 








LBQS-epoch C IV Emission EW vs. BAL EW 


332] 


56 


14 


HET-epoch C IV Emission EW vs. BAL EW 


3321 


41 


14 


Change in C IV Emission EW vs. BAL EW 


3321 


66 


14 
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Table 3. C IV BAL Absorption Equivalent Widths' 



LBQS Name 


LBQS 


Palomar 




SDSS 


HEX 


0009 + 0219 


-12.36 ± 2.07 








-30.23 ±0.76 


0019 + 0107 


-21.90 ± 1.23 


-26.35 ±0, 


.88 




-33.74 ±0.84 


0021 - 0213 


-42.91 ± 1.54 








-47.78 ±0.79 


0022 + 0150 


-15.98 ± 1.68 








-10.36 ±0.84 


0025 - 0151 


-26.99 ± 1.41 








-16.00 ±0.91 


0029 + 0017 


-31.99 ± 1.54 






-33.82 ±0.71 


-38.02 ±0.89 


1231 + 1320 


-27.34 ± 1.53 








-33.56 ±0.79 


1235 + 0857 


-26.53 ± 1.42 


-28.60 ± 0, 


.79 




-23.57 ±0.88 


1235 + 1453 


-35.64 ±3.08 


-21.86 ±0, 


.86 


-39.80 ±0.87 


-37.94 ±0.97 


1240 + 1607 


-28.14 ± 1.48 


-36.30 ± 0, 


.94 




-34.40 ± 0.79 


1243 + 0121 


-41.77 ± 2.53 


-40.38 ± 0, 


.79 


-33.63 ±0.70 


-36.85 ±0.78 


1314 + 0116 


-22.14 ±2.37 


-22.08 ±0, 


.90 


-30.22 ±0.73 


-42.10 ±0.73 


1442 - 0011 


-37.21 ± 1.24 


-33.48 ±0, 


.83 




-39.76 ±0.77 


1443 + 0141 


-43.99 ± 1.56 


-48.83 ±0, 


.81 




-44.86 ± 0.76 



^Equivalent widths are given in Angstroms. 



Table 4. Si IV And C IV BAL Absorption Equivalent Widths (Out To -20, 000 km s"^) 



LBQS Name 


LBQS Si IV 


LBQS C IV 


HET Si IV 


HET C IV 


0009+0219 


-0.79 + 1.28 


-7.74 ± 1.83 


-6.45 ±0.69 


-23.67 ±0.61 


0022+0150 


-3.58 ±1.14 


-11.01 ± 1.44 


-3.05 ±0.68 


-6.09 ±0.70 


1231+1320 


-29.32 ± 1.04 


-10.31 ± 1.32 


-21.05 ±0.66 


-14.64 ±0.70 


1235+0857 


-5.15 ± 1.00 


-22.39 ± 1.15 


-5.55 ±0.74 


-21.20 ±0.68 


1235+1453 


-1.08 ±1.95 


-29.63 ± 2.73 


-8.23 ±0.87 


-29.55 ±0.79 


1240+1607 


-11.53 ± 1.08 


-25.40 ± 1.20 


-9.96 ±0.73 


-30.24 ±0.59 


1243+0121 


-8.90 ±1.59 


-39.92 ± 2.11 


-12.12 ±0.66 


-32.24 ±0.60 


1314+0116 


-3.81 ± 1.36 


-21.97 ±2.06 


-9.55 ±0.68 


-34.93 ±0.56 


1443+0141 


-12.00 ± 1.19 


-32.29 ± 1.27 


-12.46 ±0.74 


-33.00 ±0.61 



^BQS and HET equivalent widths measured from to —20,000 km s ^, given in 
Angstroms. 
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Table 5. C IV Emission Equivalent Widths' 



LBQS Name LBQS Palomar SDSS HET 



0009 


+ 0219 


18.23 








9.79 


0019 


+ 0107 


32.51 


26, 


.79 




29.87 


0021 


- 0213 


17.48 








16.67 


0022 


+ 0150 


31.58 








23.49 


0025 


- 0151 


34.11 








22.36 


0029 


+ 0017 


58.59 






44.79 


50.35 


1231 


+ 1320 


22.76 








8.92 


1235 


+ 0857 


55.72 


25, 


.92 




29.89 


1235 


+ 1453 


22.35 


33 


.56 


45.01 


35.25 


1240 


+ 1607 


29.53 


37, 


.00 




19.03 


1243 


+ 0121 


44.80 


39, 


.46 


45.44 


46.15 


1314 


+ 0116 


57.99 


53, 


.78 


41.11 


31.35 


1442 


- 0011 


20.88 


28, 


.27 




25.00 


1443 


+ 0141 


22.77 


29, 


.92 




21.19 



''AH EWs are given in Angstroms. We do not include 
formal errors because the uncertainty is primarily a fac- 
tor of our manual adjustments to generate consistent 
fits for a variety of complex and absorbed line profiles 

(E3. 
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Fig. 1. — LBQS, Palomar, SDSS, and HET spectra (when available) plotted as flux density 
(erg s~^ cm~^ A~^), with an arbitrary y-axis normalization for each panel. The source name 
is shown at top (e.g., "0009+0219") and observing epochs are listed for each panel. Spectral 
properties for each observing epoch and model fitting techniques are described in §2.11 The 
spectra have been smoothed with a boxcar filter 3 bins wide for visual clarity. The thick red 
line indicates our emission model fit in each case. Error bars have been omitted for clarity, 
but are representative of the scatter in each spectrum. A color version and an extended 
version of this figure are available in the online journal. 
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< 

c 
Q 



X - 



Palomar (black), LBQS (red) j 




Palomar (black), SDSS (red) 




Palomar (black), HET (red) J 




1500 2000 
Wavelength (A) 



2500 



Fig. 2. — Comparison of the spectra obtained for 1235+1453. The spectra have been 
smoothed using a boxcar of width 7 bins for visual clarity. In each panel, the Palomar 
epoch is shown in black, while the red line represents the LBQS (top), SDSS (middle) and 
HET (bottom) spectra. The shaded portion represents the potential C IV BAL absorption 
region. A color version of this figure is available in the online journal. 
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00 
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■ I—' 
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- 
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1000 
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Fig. 3. — The Palomar spectrum of 1235+1453 (binned as in Figure |2]), with two alternative 
models for continuum, Si IV, and C IV emission plotted as thick lines. The top model is the 
one we use in this work, selected to match BAL absorption in other epochs more closely. 
The bottom model attempts to match the continuum underlying the emission line at 1300 A, 
but falls below the emission in the 1400-1450 A range and in the absorbed region blueward 
of 1200 A. See ^2.21 for further discussion. 
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1300 
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1450 



Rest-Frame Wavelength (A) 



Fig. 4. — The Si IV BAL absorption regions of LBQS (black), Palomar (red, when available), 
SDSS (green, when available), and HET (blue) spectra for our sources, with continuum fits 
divided out. The spectra have been smoothed by a boxcar filter 3 bins wide for visual clarity. 
Typical error bars for the unsmoothed spectra are shown at the bottom right, color-coded 
for each epoch. A color version and an extended version of this figure are available in the 
online journal. 
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Fig. 5. — Same as Figure HI but for the C IV absorption region. Thick horizontal hues (color- 
coded to match the spectra) indicate regions that are formally identified as components of 
a C IV BAL in the calculation of BIq (defined in §3.2p . Because the formal definition of 
the balnicity index only extends out to outflow velocities —25,000 km s~^, no absorption 
blueward of ^1A2Q A is included. A color version and an extended version of this figure are 
available in the online journal. 
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1000 2000 

Rest- Frame At (days) 

Fig. 6. — The change in EW over time for sources in our study. Because absorption EW is 
negative, strengthening BALs are in the lower half of the plot. Here, EWs are calculated 
by integrating over all bins in the velocity range -30,000 to km s"^ from the C IV rest 
wavelength. The observing date of the first (LBQS) epoch is defined to be Atsys = for each 
source, and the EW of the first (LBQS) epoch is subtracted from all EWs so that only the 
change in EW is shown. Lines connect measurements corresponding to a single source, and 
each individual measurement is marked with a letter indicating the corresponding observing 
campaign ("P" for Palomar, "S" for SDSS, and "H" for HET). The square root, s, of the 
unbiased sample variance (in units of 1000 km s~^) is shown for comparison between sources 
with Atsys < 500 and Atsys > 2000 days. In the color version (online), the points and lines 
for each source have been given different (arbitrary) colors for visual clarity. A color version 
of this figure is available in the online journal. 
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Fig. 7. — Histogram (thin solid line) of AEW over 3-7 yr for the combined sample of 23 
sources from this work and those of G08 with duplicates removed. We calculate AEW 
as the difference in EW between the earliest (LBQS) and latest (SDSS or HET) epochs. 
Strengthening BALs have AEW < 0. For comparison, we also plot a histogram (dotted 
line) of AEW over time scales of weeks to months between the LBQS and Palomar epochs 
in our sample, when available. 
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Fig. 8. — Magnitude of fractional change in EW against the rest-frame time between epochs. 
Points marked with "P" correspond to Palomar epochs, when available, while points marked 
with "H" correspond to HET epochs. Solid lines connect Palomar and HET epochs for the 
same source. In the color version (online), the points and lines for each source have been 
given different (arbitrary) colors for visual clarity. A color version of this figure is available 
in the online journal. 
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Fig. 9. — Comparison of rest-frame AEW/At measured for the LBQS-to-Palomar epochs 
(x-axis) and the Palomar-to-HET epochs (y-axis). The time duration between the Palomar 
and HET epochs is much longer, leading to a much smaller range of values on the y-axis and 
larger error bars in the x-direction. The dotted line indicates the line x = y. 



-se- 



ns 

DC 



in 



LO 

o 



^ o 
O 

< o 

O LO 

> ? 



\ Si' IV 0069+0219' 




LBQS ; 

HET - 


— 1 — 1 — 1 — 1 — — 1 — 1 — 1 — 1 — — 1 — 1 — 1 — 1 — — h 

; C IV 


— 1 — 1 — 1-^ 




TBQS^ 
HET - 


H — 1 — 1 — 1 — — 1 — 1 — 1 — 1 — — 1 — 1 — 1 — 1 — — 1- 
- 


r' — ' ' - 

Si \V r 

C IV ; 




x104 -1.5x104 -104 -5000 
Outflow Velocity (1000 km/s) 



^-^x104 -1.5x104 -104 -5000 
Outflow Velocity (1000 km/s) 



Fig. 10. — For our two candidate cases for coordinated Si IV and C IV absorption variation 
(0009+0219, on left, and 1443+0141, on right), we show the ratio spectra (with emission 
model divided out) for Si IV (top panel) and C IV (middle panel) absorption regions out to 
—20000 km s"^. LBQS epochs are shown with thin lines and HET epochs are shown with 
thick lines in those two panels. In the bottom panel, we plot AS{v) as a thin line and AC{v) 
as a thick line. These "difference spectra," described in §3.5[ represent the difference between 
the absorber transmission for Si IV and C IV, respectively, as a function of velocity. Related 
changes in the absorber strength can be seen in at least a subset of velocity regions for each 
ion (e.g., increased absorption at -8000 and -13,000 km s~^ and decreased absorption at 
-4000 km s-i for 1443+0141). 
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Fig. 11.— Change in EW (AEW) between LBQS and HET epochs for Si IV and C IV 
absorption regions for sources at sufficiently high redshift to have full HET coverage of the 
Si IV absorption region. The Si IV and C IV EWs are integrated out to —20, 000 km s~^. 




|A EW / <EW>| for C IV 

Fig. 12. — Same as Figure [11] but showing the fractional change in EW for Si IV (y-axis) 
and C IV (x-axis). The sohd hne represents a fit to the data with the outher 1235+1453 
removed. 



